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Abstract. We connect quantum graphs with infinite leads, and turn them to scat- 
tering systems. We show that they display all the features which characterize quantum 
scattering systems with an underlying classical chaotic dynamics: typical poles, delay 
time and conductance distributions, Ericson fluctuations, and when considered sta- 
tistically, the ensemble of scattering matrices reproduce quite well the predictions of 
appropriately defined Random Matrix ensembles. The underlying classical dynamics 
can be defined, and it provides important parameters which are needed for the quan- 
tum theory. In particular, we derive exact expressions for the scattering matrix, and an 
exact trace formula for the density of resonances, in terms of classical orbits, analogous 
to the semiclassical theory of chaotic scattering. We use this in order to investigate the 
origin of the connection between Random Matrix Theory and the underlying classical 
chaotic dynamics. Being an exact theory, and due to its relative simplicity, it offers 
new insights into this problem which is at the fore-front of the research in chaotic 
scattering and related fields. 
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1. Introduction 



Quantum graphs of one-dimensional wires connected at nodes were introduced already 
more than half a century ago to model physical systems. Depending on the envisaged 
application the precise formulation of the models can be quite diverse and ranges from 
solid-state applications to mathematical physics Jj], p), |], U) IE H 0> H H> H0> EL El 



Lately, quantum graphs attracted also the interest of the quantum chaos community 
because they can be viewed as typical and yet relatively simple examples for the large 
class of systems in which classically chaotic dynamics implies universal correlations in 
the semi classical limit |U| [14|, [15], [U| [17|, [18], [L9|, |2(J, ^TJ . Up to now we have only a limited 
understanding of the reasons for this universality, and quantum graph models provide a 
valuable opportunity for mathematically rigorous investigations of the phenomenon. In 
particular, for quantum graphs an exact trace formula exists [12, [13, [L4j] which is based 



on the periodic orbits of a mixing classical dynamical system. Moreover, it is possible 
to express two-point spectral correlation functions in terms of purely combinatorial 
problems |16|, \V7\, [18], |19| . 

By attaching infinite leads at the vertices, we get non-compact graphs, for which a 
scattering theory can be constructed [^, ^3], ^4|. They display many of the feature 
that characterize scattering systems with an underlying chaotic classical dynamics 
25| , |26| , [2T| , [28"| , p9[ , and they are the subject of the present paper. The quantum 
scattering matrix for such problems can be written explicitly, together with a trace 
formula for the density of its resonances. These expressions are the analogues of the 
corresponding semiclassical approximations available in the theory of chaotic scattering 
]], albeit here they are exact. With these tools we analyze the distribution 



26 



of resonances, partial delay times, and the statistics of the fluctuating scattering 
amplitudes and cross sections. Moreover, we address issues like the statistical properties 
of the ensemble of the scattering matrices, and conductance distribution. Finally 
we analyze the effect of non-uniform connectivity on the statistical properties of the 
scattering matrix. A main part of our analysis will be focused on the comparison of the 
statistical properties of the above quantities with the Random Matrix Theory (RMT) 
predictions pi, [31], pi pi M, M, M, E7L gBL [391 |0 . 



The paper is structured in the following way. In section (|2]), the mathematical model 
is introduced and the main definitions are given. Section (|3]) is devoted to the derivation 
of the scattering matrix for graphs. The trace formula for the density of resonances 
of quantum graphs is presented in section (|4]) which includes also the analysis of the 
underlying classical system. The next section, is dedicated to the analysis of various 
statistical properties of the 5*-matrix. Our numerical data are compared both with the 
predictions of RMT and with the semiclassical expectations. Finally, in section (^), we 
analyze the class of star-graphs for which several results can be analytically derived. 
Our conclusions are summarized in section (17|). 
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2. Quantum Graphs: Definitions 

We start by considering a compact graph Q. It consists of V vertices connected by B 
bonds. The number of bonds which emanate from each vertex i defines the valency Vi 
of the corresponding vertex (for simplicity we will allow only a single bond between 
any two vertices). The graph is called v — regular if all the vertices have the same 
valency v. The total number of bonds is B — | J2i=i v i- Associated to every graph is its 
connectivity matrix C . It is a square matrix of size V whose matrix elements Cij take 
the values 1 if the vertices i,j are connected with a bond, or otherwise. The bond 
connecting the vertices % and j is denoted by b = and we use the convention that 

% < j. It will be sometimes convenient to use the "time reversed" notation, where the 
first index is the larger, and b = (j, i) with j > i. We shall also use the directed bonds 
representation, in which b and b are distinguished as two directed bonds conjugated by 
time-reversal. We associate the natural metric to the bonds, so that x it j (xj ti ) measures 
the distance from the vertex % (j) along the bond. The length of the bonds are denoted 
by L b and we shall henceforth assume that they are rationally independent. The mean 
length is defined by (L) = (1/-B) Y2b=i ^ b ano - ^ n an numerical calculations bellow it will 
be taken to be 1. In the directed- bond notation L b = L~ h . 

The scattering graph Q is obtained by adding leads which extend from M(< V) 
vertices to infinity. For simplicity we connect at most one lead to any vertex. The 
valency of these vertices increases to — Vi + 1. The M leads are denoted by the index 
i of the vertex to which they are attached while Xi now measures the distance from the 
vertex along the lead i. 

The Schrodinger operator (with h~ = 2m = 1) is defined on the graph Q in the 
following way: On the bonds b, the components of the total wave function \l> are 
solutions of the one - dimensional equation 

{'^~ Ab ) ^b(x) = k 2 ^ b (x), b=(i,j) (1) 

where A b (with lR.e(A b ) ^ and A b = —A^) is a "magnetic vector potential" which breaks 
the time reversal symmetry. In most applications we shall assume that all the A b s are 
equal and the bond index will be dropped. The components of the wave functions on 
the leads, ^i(x), are solutions of 

- ^2*i0*0 = k 2 ^i{x), i = l,..,M. (2) 

At the vertices, the wavefunction satisfies boundary conditions which ensure current 
conservation. To implement the boundary conditions, the components of the wave 
function on each of the bonds b and the leads i are expressed in terms of counter 
propagating waves with a wave-vector k: 

On the bonds : ^ b = a b e i{k+Ab)xb + c b e t{ ~ k+Ab)xb 

On the leads : % = ^e"^ + . (3) 
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The amplitudes a b , c b on the bonds and I i} on the lead are related by 
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These equalities impose the boundary conditions at the vertices. The vertex scattering 
matrices are Vi x vi unitary symmetric matrices, and j,j' go over all the V{ bonds 
and the lead which emanate from i. The unitarity of £W guarantees current conservation 
at each vertex. 

On the right hand side of (^), the vertex scattering matrix £w was written 
explicitly in terms of the vertex reflection amplitude the lead-bond transmission 



amplitudes {t^}, and the Vi x vi bond-bond transition matrix erf'-,, which is sub unitary 



3,3 



(| detail < 1), due to the coupling to the leads. 

Graphs for which there are no further requirements on the £W shall be referred 
to as generic. It is often convenient to compute the vertex scattering matrix from 
a requirement that the wave function is continuous and satisfies Neumann boundary 
conditions at that vertex. These graphs shall be referred to as Neumann graphs, and 
the resulting £W matrices read: 

2 (j) 2 2 



~(i) 



(i) 



(5) 



Vertices which are not coupled to leads have p l = 1, r„- = 0, while the bond-bond 



transition matrix a^j, 



— 5a ji 

v 3,3 



is unitary. 



3. The S-matrix for Quantum Graphs 

It is convenient to discuss first graphs with leads connected to all the vertices M = 
V. The generalization to an arbitrarily number M < V of leads (channels) is 
straightforward and will be presented at the end of this section. 

To derive the scattering matrix, we first write the bond wave functions using the 
two representations which are conjugated by "time reversal": 

^ b (x b ) = a b e i{k+Ab)xb + c b e i{ ~ k+Ab)xb = 
\& b (x h ) = a b e i(k+A i> )x -» + c b e i{ - k+A ^ = 

= aye i ^~ k ~ A ' b)Xb e i ^ k+A ' b)Lb + cys i ^~ k+A ' b)Lb e i ^ k ~ A ' b)Xb . (6) 

Hence, 

c b = a b e^ k+A ~ ALb , a b = c b e^- k+A ^ Lb . (7) 

In other words, but for a phase factor, the outgoing wave from the vertex % in the 
direction j is identical to the incoming wave at j coming from i. 



Quantum Graphs: A simple model for Chaotic Scattering 



5 



Substituting a b from Eq. (0) in Eq. (f|), and solving for Cjj we get 

= ^)/ i + ^ r j ; %, (8) 
f 

where 1 is the 2B x 2B unit matrix. Here, the "bond scattering matrix" Sf, is a sub- 
unitary matrix in the 2B dimensional space of directed bonds which propagates the 
wavef unctions. It is defined as Ssik, A) = D(k;A)R, with 

D ijti/f (k, A) = S i>v S jtjl e ikL ^^ (9) 

^jifTim (J n,i x ~ J j,i'^ / i,mVji t i m ' 

D(k,A) is a diagonal unitary matrix which depends only on the metric properties of 
the graph, and provides a phase which is due to free propagation on the bonds. The 
sub-unitary matrix R depends on the connectivity and on the bond-bond transition 
matrices a. It assigns a scattering amplitude for transitions between connected directed 
bonds. R is sub-unitary, since 

v 

|det£| = JJ | deta (i) | < 1. (10) 

i=l 

Replacing c it ji in the second of Eqs. (H) we get the following relation between the 
outgoing and incoming amplitudes Oi and Ij on the leads: 

O^f^h + Y^Tf (l-S B (k;A)Y\ D^r^Ii (11) 

^ \ / (i,r),(s,i) 



3 jrs 



Combining (|TT|) for all leads % = 1, . . . , V, we obtain the unitary V x V scattering matrix 

S (v )} 

S% ] = S itjP ® + £ r« (l - S B (*; A)) ~\ D (sj) t^. (12) 

r,s 

From Eq. (P), we see that the scattering matrix may be decomposed into two parts 
S( v \k) = S D + S-f l (k) which are associated with two well separated time scales of the 
scattering process. S D (= o~i,jP l ) is the prompt reflection at the entrance vertex and 
induces a "direct" component. In general, it varies very slowly with energy, and is 
system dependent. On the other hand the "chaotic" component of the S matrix, S^ l (k), 
starts by a transmission from the incoming lead % to the bonds (i,r) with transmission 
amplitudes Tr . The wave gains a phase e l ( k+A b) L b f or eac h bond it traverses and a 
scattering amplitude aifs at each vertex. This multiple scattering inside the interaction 
region becomes apparent when the expansion 

oo 

(l-Sni^A^-^J^^A) (13) 

n=0 
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is substituted in (|i~2p. Eventually the wave is transmitted from the bond {s,j) to the 
lead j with an amplitude Ts • Explicitly, 

S$P = Stjp® + Bte i{kk+et) (14) 

where is the set of the trajectories on Q which lead from i to j. Bt is the 

amplitude corresponding to a path t whose length and directed length are l t = J2b&t 
and & t = ^2b€t LbAb respectively. Thus the scattering amplitude S^p is a sum of a 
large number of partial amplitudes, whose complex interference brings about the typical 
irregular fluctuations of |<S^| 2 as a function of k. 

One of the basic concepts in the quantum theory of scattering are the resonances. 
They represent long-lived intermediate states to which bound states of a closed system 
are converted due to coupling to continua. On a formal level, resonances show up as 
poles of the scattering matrix occurring at complex wave-numbers K n = k n — |T n , 
where k n and T n are the position and the width of the resonances, respectively. From 
(|j~2|) it follows that the resonances are the complex zeros of 

(c(K) = det(l-S B (K;A)) = 0. (15) 



The eigenvalues of Sb are in the unit circle, and therefore the resonances appear in 
the lower half of the complex k plane. Moreover from Eq. fll~5|) it is clear that their 
formation is closely related to the internal dynamics inside the scattering region which 
is governed by Sb- 

There exists an intimate link between the scattering matrix and the spectrum of the 
corresponding closed graph. It manifests the exterior -interior duality [f|l]] for graphs. 



The spectrum of the closed graph is the set of wave-numbers for which S™' has +1 as 
an eigenvalue. This corresponds to a solution where no currents flow in the leads so 
that the conservation of current is satisfied on the internal bonds. 1 is in the spectrum 
of S (v ^ if 

( g {k) =det [l-S iv) {k)] =0 . (16) 
Eq. (H) can be transformed in an alternative form 



(r) (r) 

Cc(fc) =det[l-p] det ^ 1 ~ jP(A:) ^ =0 ; R ir . si = Ri r . sj + 5 rs - ^ , . (17) 

which is satisfied once 

det[l - D(k)R] = 0. (18) 

In contrast to R, R is a unitary matrix in the space of directed bonds, and therefore 
the spectrum is real. flTSP is the secular equation for the spectrum of the compact part 
of the graph, and it was derived in a different way in fl3| . 



The difference 8R = R — R gets smaller as larger graphs are considered (for graphs 
with Neumann boundary conditions it is easy to see that the difference is of order -). 
That is, the leads are weekly coupled to the compact part of the graph, and one can use 
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Figure 1. Poles of the -matrix for regular Neumann graphs. The exact evaluated 
poles are indicated with (o) while (*) are the results of the perturbation theory (p0|): 
(a) V = 5 and v — 4 and (b) V = 15 and v = 14. 



perturbation theory for the computation of the resonance parameters. To lowest order, 
(SR = 0), the resonances coincide with the spectrum of the compact graph. Let k n be 
in the spectrum. Hence, there exists a vector \n) which satisfies the equation 

D{k n )R \n) = 1 \n) . (19) 

To first order in 5R, the resonances acquire a width 



5K n = —i 



(n 


D(k n )5R\n) 


(n 


L 


n) 



(20) 



To check the usefulness of this result, we searched numerically for the true poles for a 
few scattering graphs and compared them with the approximation (|20|). In figure 1 we 
show the comparison for fully connected Neumann graphs with V = 5, 15 and A = 0. As 
expected the agreement between the exact poles and the perturbative results improves 
as v increases. 

We finally comment that the formalism above can be easily modified for graphs 
where not all the vertices are attached to leads. If the vertex / is not attached, one has 
to set = 1, r. = in the definition of E^. The dimension of the scattering matrix 
is then changed accordingly. 

For the sake of completeness we quote here an alternative expression for the S 
matrix which applies for Neumann graphs, exclusively |TJ]]. For this purpose we define 
the V x V matrix 

{- Ern^i C i,m COt(kL ijm ), i = j 
(21) 
C id e- iA ^{sm{kL itj ))-\ i / ., 

in terms of which, 

S {v) = {il + h{k))-\il - h(k)) (22) 
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where 1 is the V x V unit matrix. S( v > is unitary since h{k) is hermitian. 

In the case of graphs connected to leads at an arbitrary set of M < V vertices with 
indices {ii} , I = 1, • • • , M, the M x M scattering matrix S"( M ) has to be modified in 
the following way 

= 2iW (h(k) + iW T W) _1 W T - 1 . (23) 

Here W^j = 5i t j is the M x V leads - vertices coupling matrix, ( W = 1 when M = V). 
This form of the S matrix is reminiscent of the expression which was introduced by 
Weidenmuller to generalize the Breit-Wigner theory for many channels and internal 
states. However, (P3|) is an exact expression which involves no truncations. 

It follows from (|23|) that in the present case, one can identify the poles of the S 
matrix with the zeros of 

C§(k) = det (h(K) + iW T W) = . (24) 

Its main advantage over (|i~5|) is that it involves a determinant of a matrix of much lower 
dimension. 



4. The trace formula for resonances and classical scattering on graphs 

In the spectral theory of bounded hamiltonian systems, the most fundamental object of 
study is the spectral density which consists of a sum of 5 functions at the spectral points. 
For open systems, it is replaced by the resonance density, which is defined on the real k 
line and consists of an infinite sum of Lorentzians which are centered at k n = 3?eK n and 
have width T n = —2QmK n , where K n are the complex poles of the scattering matrix. In 
the present section we shall express the resonance density for scattering graphs in terms 
of periodic orbits of their compact part. This is the analogue of the trace formula for 
bounded graphs. 



4-1. The Trace Formula 

The resonance density dn(k) can be deduced from the total phase 
4 In det S( M \k) of the scattering matrix 



d ^^§ = -^l lndets<M, w- < 25) 

It is a smooth function for real k, and can be interpreted as the mean length of the 
delay associated with the scattering at wave- number k [^q| . 

Using Eq. (|i~2| ) and performing standard manipulations we obtain the following 
expression for the phase $(fc) 

- $(0) = -29fmlndet(J- S B (k; A)) + Ck . (26) 

where £ = 2 J2b=i Lb is twice the total length of the bonds of Q. 
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Using the expansion 

oo 1 

lndet( I - S B (k; A)) = - V -trSl(k; A) (27) 

z — ' n 

n=l 

we rewrite Eq. ( p6|) as 

oo 1 

$(jfc) = $(0) + £A; + 23m V -trSSffc; A). (28) 

n=l 

On the other hand, using Eq. (|9]) we can write tnSjg(/s; A) as sums over n— periodic 
orbits on the graph 

ti (S^(h, A)) = n pA r p e i{lpk+@p)r , (29) 
peVn 

where the sum is over the set V n of primitive periodic orbits whose period n p is a divisor 
of n, with r = n/n p (primitive periodic orbits are those which cannot be written as 
a repetition of a shorter periodic orbit). The amplitudes A p are the products of the 
bond-bond scattering amplitudes a^\, along the primitive loops i.e. 

rip 

A = W~^l- (30) 

i=l 

Substituting Eqs. (^,^9|) in Eq. (p3[) one gets the resonance density 

oo 

d R (k) = —C + -^ e y Y n p l p rAy {l * k+e v )r (31) 

n=l pGPn 

Eq. (^) is an exact trace formula for the resonance density. The first term on the right 
hand side of Eq. ([JT]) corresponds to the smooth resonance density, while the second 
provides the fluctuating part. We notice that the mean resonance spacing is given by 
2tt 2tt 

A = — ~ (32) 

C 2B(L) K ' 

and it is the same as the mean level spacing obtained for the corresponding bounded 
graph Q 



4-2. Classical Dynamics 

We conclude this section with the discussion of the classical dynamics on the graph Q. 
A classical particle moves freely as long as it is on a bond. The vertices are singular 
points, and it is not possible to write down the analogue of the Newton equations there. 
In []T3 , 14]] it was shown that it is possible to define a classical evolution on the graph: 



A Poincare section on the graph consists of the discrete set of directed bonds. The 
phase-space density at a (topological) time n is the set of occupation probabilities p&(n) 
of the directed bonds, and the classical evolution is governed by a Markovian master 
equation. Applied to the compact part of a scattering graph it reads, 

p b {n + 1) = ^2 Ub,b'Pb'(n) (33) 
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where the transition matrix U^y is given by the corresponding quantum transition 
probability 

Notice that U does not involve any metric information on the graph. 

Due to loss of flux to the leads Yly U},y < 1, and the phase-space measure is not 
preserved, but rather, decays in time. The probability to remain on Q is 

2B 

P{n) = J> 6 (n) =^UwPv{n - 1) ~ e" r ^P(0) (35) 

6=1 6,6' 

where exp(— r c /) is the largest eigenvalue of the "leaky" evolution operator Ubv- 

For the w-regular graph with cfW given by @ the spectrum of U is restricted to 
the interior of a circle with radius given by the maximum eigenvalue V\ — (v — l)r 2 + p 2 
with corresponding eigenvector |1 >= (1/25) (1, 1, ■ • ■ , 1) T . Hence the decay rate 
r c ; = — In v\ for regular Neumann graphs take the simple form 

r d = -ln(l-r 2 ) ^{2/{l+v)f. (36) 

We notice that removing the leads from the vertices and turning Q into a compact graph 
Q we get r c z = since in this case r^' = (and p^ = 1) and the phase-space measure 
is preserved as expected. 

The inverse decay rate T c i = T^ 1 , gives the average classical delay time that the 
particle spends within the interaction region. Injecting a particle from the leads to the 
scattering domain, its probability to be on any bond randomizes, because at each vertex 
a Markovian choice of one out of v directions is made. The longer a particle remains 
within the interaction regime, the more scattering events it experiences. The set of 
trapped trajectories whose occupancy decays exponentially in time is the analogue of 
the strange repeller in generic Hamiltonian systems displaying "chaotic scattering" . 

5. Statistical Analysis of the S— matrix 

So far we developed the scattering theory of graphs, pointing out their similarity with 
scattering systems which display chaotic scattering in the classical limit. Due to the 
interference of a large number of amplitudes, the S— matrix fluctuates as a function of 
k, and its further analysis calls for a statistical approach which will be the subject of 
the present section. We shall show that quantum graphs possess typical poles, delay 
time and conductance distributions, Ericson fluctuations of the scattering amplitudes, 
and when considered statistically, the ensemble of scattering matrices are very well 
reproduced by the predictions of RMT. At the same time deviations from the universal 
RMT results, which are related to the system- specific properties of some graphs, will 
be pointed out. The study of these deviations is especially convenient for graphs 
because of the transparent and simple scattering theory developed in terms of scattering 
trajectories. 
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Figure 2. Representative examples of scattering cross sections |5£p| 2 for Neumann 
graphs in a regime of (a) isolated resonances (V — 15, v = 14 with corresponding 
(7), = 0.59) and (b) overlapping resonances (V = 49, v = 14 with corresponding 
(7). ~ 2). The real parts of the resonances in this energy interval are indicated by the 
vertical lines. 

An important parameter which is associated with the statistical properties of the 
S-matrix is the Ericson parameter defined through the scaled mean resonance width as: 

<7> fc - (37) 

where denotes spectral averaging and A is the mean spacing between resonances. 
The Ericson parameter determines whether the resonances overlap ((7)* > 1) or are 
isolated ((7)* < 1). Typical example for the two extreme situations are shown in 
figure 2. 

The degree of resonance overlap determines the statistical properties of the S- 
matrix. We shall show in the sequel that the mean width can be approximated by the 
classical decay rate (7)^ = 7^. For the v regular graphs discussed above, we have 

T d 4 v V 
Id = -T- ~ 7T TT n • 38 

A Z7T 1 + vl + V 



where we made use of Eqs. ( ]32]J36| ). Thus changing v and V we can control the degree 
of overlap allowing to test various phenomena. 

In what follows, unless explicitly specified, we shall consider regular graphs with 
one lead attached to each vertex, i.e. M — V. Finally, the widths are always scaled by 
the mean spacing A i.e. 7 n = 
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5.1. The resonance width distribution 

The resonance width distribution can be computed for a given graph in the following 
way. Consider the complex k = x + iy plane, where the zeros of the secular function 

Ce(«)=det(l-5 B («;A)) = f(x,y) + ig(x,y) . (39) 

are the poles of the S matrix (resonances). The variables (x, y) are expressed in units of 
A. On average, the number of resonances with real part, x G [0, X] is X. The density 
of resonance widths reads 



dn 



dx . (40) 



1 f A 

V(i) = lim — / 5 (f(x, y = -7)) 5 (g(x, y = -7)) 

We used the Cauchy-Riemann theorem for the evaluation of the Jacobian (f x g y — f y g x ) 
which multiplies the two 5 functions that locate the complex zeros of (g. We now recall 
that the k dependence of Sb{k~, A) comes from the factors e tKLb in (0). This implies that 
for a given y, (g is a quasi periodic function of x. Moreover, expanding the determinant 



(39) it is not difficult to show that the frequencies involved can be written as linear 
combinations of the bond lengths A = J^QbLb with integer coefficients g& = 0,1, or 2. 
Since the bond lengths are rationally independent, we find that (g depends on a finite 
number of incommensurate frequencies. The expression fl4*0"D can be regarded as a "time" 
integral over a trajectory on a multidimensional incommensurate torus, which covers the 
torus ergodically. Hence, the integral can be replaced by a phase space average, 

2 



dn 



(41) 



Here ip denotes the vector of independent angles on the J dimensional torus. Although 
the above formula provides a general framework, its application to actual graphs is a 
formidable task. 

An important feature of the distribution of the resonances in the complex plane 
can be deduced by studying the secular function (g(n). Consider Cg{ K = 0). If one of 
the eigenvalues of the matrix R (g) takes the value 1, (g(k = 0) = and because of the 
quasi-periodicity of (g, its zeros reach any vicinity of the real axis infinitely many times. 
The largest eigenvalue of the R matrix for w-regular Neumann graphs is 1 , and therefore 
the distribution of resonance widths is finite in the vicinity of 7 = 0. For generic graphs, 
the spectrum of R is inside a circle of radius X m ax < 1- This implies that the poles are 
excluded from a strip just under the real axis, whose width can be estimated by 

Fgap — — 21n(|A maa . \)/L max . (42) 

where L max is the maximum bond length. The existence of a gap is an important feature 
of the resonance width distribution V(j) for chaotic scattering systems. 



A similar argument was used recently in [23] in order to obtain an upper bound for 
the resonance widths. It is 

Tmax ~ — 2 ln(|A m j n |)/L TO j n , (43) 
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Figure 3. (a) The 5000 resonances of a single realization of a complete V = 5 graph 
with A ^ and M — V. The solid line marks the position of the gap 7 sa p- (b) The 
distribution of resonance widths V^y). The solid line is the RMT prediction (EH). The 
difference 7^(7) — Vcue{i) is shown in the inset. 



where X min and L min are the minimum eigenvalue and bond length, respectively. 

The distribution of the complex poles for a generic fully connected graph with V = 5 
is shown in figure 3a. The vertical line which marks the region from which resonances 
are excluded was computed using (|4~^). 

Random matrix theory can provide an expression for the distribution of resonances. 
In the case of non-overlapping resonances, perturbation theory shows that the resonance 
widths are distributed according to the so-called x 2 distribution 

ni) = ? M lTo^L (A-) PM2 1 exp(- 7 /3M/2( 7 ) fc ), (44) 



( 7 > fcJ T(/W/2) V<7>/ 

where /3 — 1(2) for systems which respect (break) time-reversal symmetry, and r(x) is 
the gamma-function. Once (7)^ becomes large enough the resonances start to overlap, 
and (f0|) does not hold. In the general case, Fyodorov and Sommers |36|, 57[] proved that 
the distribution of scaled resonance widths for the unitary random matrix ensemble, is 
given by 

l7j r(M) 7 d^ M V M ) 

where the parameter g = ^z^m^\ ~ 1 controls the degree of coupling with the channels 
(and is assumed that ^ = g Wi = 1,...M). For g ^> 1 (i.e. weak coupling regime) 
Eq. (|5D reduces to (||). 

In the limit of M > 1, Eq. (^) reduces to the following expression p7 

otherwise 
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Figure 4. Resonance width distribution V("f) for a complete Neumann graph with 
M = V = 5 and A ^ 0. The solid line is the RMT prediction (g§). 

It shows that in the limit of large number of channels there exist a strip in the complex 
k— plane which is free of resonances. This is in agreement with previous findings 
p7| , |30| , fPfl . In the case of maximal coupling i.e. g — 1, the power law (f4~6|) extends to 
infinity, leading to divergencies of the various moments of 7's. Using fl4"5| ) we recover 
the well known Moldauer-Simonius relation |Kj for the mean resonance width [[37] 



< 7>l = _2^M. (47) 

The resonance width distribution for a V = 5 regular and generic graph is 
shown in figure 3b together with the RMT prediction, which reproduces the numerical 
distribution quite well. Figure 4 shows a similar comparison for a Neumann graph. The 
relatively high abundance of resonances in the vicinity of the real axis conforms with 
the expectations. 



5.2. The Form Factor 

To investigate further the dynamical origin of the resonance fluctuations, we study the 
resonance two-point form factor K(t). The main advantage of K(t) is that it allows us 
to study the resonance fluctuations in terms of classical orbits. It is defined as 

K R (t) =jd X e l2 ^ a R 2 ( X ). (48) 

where t measures lengths in units of the Heisenberg length Ijj = C and Rq,{x) is the 
excess probability density of finding two resonances at a distance Xi 

R 2 (x; ko) = (d R (k + |) d R (k - = (A d R {k + |) d R (k - ^)dk) kQ .(49) 
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Figure 5. The form factors Kn{t) for a complete V = 5 graph, with either generic 
or Neumann boundary conditions, 4^0 and M = V . The data were averaged over 
5000 spectral intervals and smoothed on small t intervals. Upper inset: Ku(t) for 
small times. Solid line correspond to the numerical data for the pentagon with generic 
boundary conditions while dashed line is the approximant Kn(t) ~ K (t)exp(—j c lt) ■ 
Lower inset: Kr{£) calculated with high resolution. The lengths of periodic orbits of 
the close graph are indicated by arrows. 



Above (■ ■ -)k indicate averaging oven a number of spectral intervals of size = 2ko, 
centered around k . Here d R {k) is the oscillatory part of d R (k) (see Eq. (|31"|)). 
Substituting the latter in Eqs. (f48l , f49|) we obtain Knit) in terms of periodic orbits and 
their repetitions 

2 



K R {t) 



2M 



(50) 



where 5 N {x) = {sin(N x / 2)) / (N x / 2) and N = ^. 
the spectral form factor of the compact graph [13|, |14 



A similar sum contributes to 
. However, the corresponding 
amplitudes are different due to the fact that A p contains also the information about the 
escape of flux to the leads. 

Assuming that all periodic orbits decay at the same rate, one would substitute 
A p with A p exp(—n p r y c i/2) where A p is the weight assign to the periodic orbit of the 
corresponding close system. Then Eq. ( |5*0"D takes the following simple form 

K R (t) w K(t)e-^ 1 , (51) 

where K{t) is the form factor of the compact system ||29|| . Notice that for t 7^ 
the resonance form factor K R {t) is equal to K(t). This is reasonable since an open 
system cannot be distinguished from a closed one during short times. This simple 
approximation is checked in the inset of figure 5 (see dashed line) and it is shown to 
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reproduce the numerical data rather well in the domain t < 5. The asymptotic decay is 
dominated by the resonances which are nearest to the gap, and it cannot be captured 
by the crude argument presented above. For generic graph, Kr decays exponentially 
but with a rate given by ^ as = ■jgap (the best fit, indicated in figure 5 by the dashed line, 
give 7 as which agrees with 7 gaj3 within 30%). For the graph with Neumann boundary 
conditions, 7 9ap = and one expects an asymptotic power-law decay. The corresponding 
best fit (see dashed line in figure 5) shows t~ 2 . 

In Hamiltonian systems in more than one dimension, the size of the spectral interval 
Afc where the spectral average is performed is limited by the requirement that the 
smooth part of the spectral density is approximately constant. Here instead we can 
take arbitrarily large spectral intervals since the smooth spectral density is constant 
|13| , |l4jl . This way, one can reach the domain where the function Kpt(t) is composed of 
arbitrarily sharp spikes which resolve completely the length spectrum for lengths which 
are both smaller and larger than In the inset of figure 5 we show the numerical K(t) 
calculated with high enough resolution. In the same inset we mark with arrows the 
location of the lengths of short periodic orbits. We notice that as long as tC is shorter 
than the length of the shortest periodic orbit, K(t) = 0. With increasing t, the periodic 
orbits become exponentially dense and therefore the peaks start to overlap, giving rise 
to a quasi continuum described approximately by Eq. (|5T|). 



5.3. Ericson fluctuations 

As was mentioned above, the scattering cross sections are dominated by either isolated 
resonances, or by overlapping resonances whose fluctuations follow a typical pattern. 
These patterns were first discussed by Ericson [|5[] in the frame of nuclear physics and 



were shown to be one of the main attributes of chaotic scattering |26| . The transition 
between the two regimes is controlled by the Ericson parameter (|37|). Typical fluctuating 
cross sections are shown in figure 2. 

A convenient measure for Ericson fluctuations is the autocorrelation function 

., Jmax 
J — Jmin 



where Aj = j max — j min + 1. To evaluate Eq. (|52|) we substitute the expression of 



the S'-matrix from Eq. (JT4]) and split the sum over trajectories into two distinct parts: 
the contributions of short trajectories are computed explicitly by following the multiple 
scattering expansion up to trajectories of length l ma x- The contribution of longer orbits 
are approximated by using the diagonal approximation, which results in a Lorentzian 
with a width -y Er . Including explicitly up to n = 3 scattering events we get, 



Jmax 

C(x; v) ~ Ge ilmaxX lEr + — V [tV xL ^+" + T 4 p 4 e 3ixL ^+- 

J — Jmin 

_|_ r 6 e *X(ij,m+im,j + ,y)j ^53^ 




Figure 6. The autocorrelation function C(x, v = 1) for regular graphs with 
Neumann boundary conditions, (o) correspond to a graph with 7 = 0.59 (isolated 
resonance regime), (*) to a graph with 7 = 1.36 (intermediate regime) while (o) to a 
graph with 7 ~ 2 (overlapping resonances regime). The solid lines correspond to the 
theoretical expression (p3|): (a) The real part of C(x, v = 1) ; (b) The imaginary part 
ofC(x^=l). 

where G is determined by the condition C(x = 0; v) = 1. The interplay between the 
contributions of long and short periodic orbits is shown in figure 6. For overlapping 
resonances, the autocorrelation function is well reproduced by a Lorentzian while for 
the case of isolated resonances one can clearly see the contributions of short paths. 

In figure 7 we report the mean resonance width (7}^ calculated numerically for 
various graphs, the parameter 7E r extracted from the best fit of the numerical C(x) 
with Eq. (|53|) ) together with the RMT prediction Eq. (|47|), and the classical expectation 
given by Eq. fl38l) . The results justify the use of the classical estimate especially in the 
limit V — > 00 for fixed v/V (which is the analogue of the semiclassical limit). In this 
limit, the RMT and the classical estimate coincide. 



5.4- S— matrix statistics 

One can check further the applicability of RMT by studying the entire S^— matrix 
distribution function. The probability density oi MxM unitary S^^-matrices is defined 
in a M + (3M(M — l)/2 parameter space and was found first in |U| to be given by the 
Poisson Kernel 

dVs(S^) = Ps(SW(S) = C ^^^^ d^S) (54) 
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Figure 7. The mean resonance width (7)^, autocorrelation width jEr, the classical 
expectation j c i, and the RMT prediction ( ft7| ) vs. for various graphs with Neumann 
boundary conditions and constant valency v — 14. 



where Cp is a normalization constant which depend on the symmetry class. All system- 
specific relevant informations are included in the ensemble average S^— matrix, 
defined as Sij = (S D )k- For regular graphs, S is proportional to the unit matrix i.e. 
Si j = fr^ij, while the eigenvectors are distributed uniformly and independent from 
the eigenphases. The invariant measure dnp{S) is given as 

M 

dfi P (S (M) ) = Yl |e* - e^f Y[ dfrdtt (55) 

i<j i=l 

where dVL is the solid angle on the M dimensional unit hyper-sphere. 

For large M values, a numerical check of this probability distribution is prohibitive. 
However, for M = 2 one can find easily the exact form of ([54]) and compare with the 
numerical results. The resulting distribution of the eigenphases 4>i,(p2 of the 2x2 
scattering matrices of regular graphs is given by 

dPs(S^) = (56) 
- ^dudip 



[1 -45cos(f)cos^ + 25 2 (2cos 2 (f) + cos(2^)) - 4S 3 cos^ cos(f ) + S 4 ] ' 



2 



where we used the notation ip = ^ 1 ^ 2 and uo = <fii — 2 - Integrating Eq. 



2 



ip we get the corresponding distribution functions V§(ui) and V§(ip). Our numerical 
results for an ensample of S'-matrices calculated for different realizations of the lengths 
of the graphs are reported in figure 8 together with the RMT predictions ( p^ ) for a 
regular graph with two channels. An overall good agreement is seen both for A = and 
A ^ 0. 




Figure 8. The V(ip) and V(w) distributions for a 2 x 2 scattering matrix S. The 
solid lines correspond to the predictions of RMT (|5^). The upper panel correspond to 
A = while the lower panel to A ^ 0. 



5.5. Partial delay times statistics 



The Wigner delay time captures the time-dependent aspects of quantum scattering. It 
can be interpreted as the typical time an almost monochromatic wave packet remains 
in the interaction region. It is defined as 

M 



2iM 



:tr 



■ S(M) f dS^(k) 
dk 



I ^-v d(j)i{k) 



2M 



dk 



(57) 



dk 



arc 



where <pi are the eigen-phases of the S^-matrix. The partial derivatives 
the partial delay times and their statistical properties were studied extensively within 
the RMT |59fl . For the one-channel case it was found that the 

probability distribution of the scaled (with the mean level spacing A) partial delay 
times Ti = i s 



2tt dk 



vf{T) 



r (|) T 2+/3/2 



2tt 



V{4>) 



(5f 



where V((p) is the Poisson Kernel Eq. (54). The general case of M-equivalent open 
channels was studied in [[35|, ^] where it was found that for broken time reversal 
symmetry the probability distribution of partial delay times is 

HT) = M lT L d{T - r s [e-^HT^Vig 2 ~ 1))] (59) 

where Io(x) stands for the modified Bessel function. 

To investigate the statistical properties of the partial delay times for our system 
we had calculated V(T) for various graphs. The resulting distributions are shown in 
figure 9 together with the RMT predictions ( |58| . |5"9"1) . An overall agreement is evident. 
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Figure 9. The distribution of the scaled partial delay times T for various graphs with 
Neumann boundary conditions. The dashed lines correspond to the RMT expectation 
(H|59|): ( a ) One channel and A = 0; (b) M = V channels and A^O. 



Deviations appear at the short time regime (i.e. short orbits), during which the "chaotic" 
component due to multiple scattering is not yet fully developed |28[ . 



5.6. Conductance Distribution 

Due to the recent experimental investigation of electronic transport through mesoscopic 
devices fHH , the study of the statistical properties of the conductance gained some 
importance. For a device connected to reservoirs by leads, the Landauer- Buttiker 
formula relates its conductance G to the transmission coefficient Tq by the expression 
G = (2e 2 IK)Tq. When each lead supports one channel, the transmission coefficient can 
be written in terms of the S'W— matrix as 

M 

T G = J2\€ ) f = 1 -\€ ) f (6°) 



where j is the input channel. 

In the absence of direct processes, the distribution of conductance V(T G ) for 
arbitrarily number of channels was worked out within RMT, and the results describe 
in a satisfactory way both the numerical calculations and the experimental data (for 
a review see 0] and references therein). However in cases where direct processes 
appear significantly, one must use the Poisson's kernel ([54]) in its full generality. This is 
exactly the case with Neumann graphs since S ^ 0. The probability distribution of the 
transmission Tq is then 



V(T a ) = / S(T G 



i -J 



(61) 




Figure 10. Conductance distribution V(Tc) for a graph with M = 2 channels: 
(a) Time reversal symmetry is preserved (A = 0). The solid line is the RMT results 
Eq. ( |6l| ) where we had used Eq. ((H^) for the Poisson Kernel while the dashed line is 
the approximate expression (p2[); (b) Broken time reversal symmetry (A ^ 0). The 
solid line is the RMT result Eq. @. 



For the case M = 2 and for a diagonal S— matrix (only direct reflection) with 
equivalent channels Sn = S 2 2 = S, Eq. ( |6lD can be worked out analytically [0 in the 
limit of strong reflection and Tq -C 1. The resulting expression is: 

-1/2 



V/3=i(Tg) 



71 



2 <\-s 2 ) Tg 



T G <^(l-S 



2\2 



For 



2 one can compute in a close form the whole distribution [|(| 

<l - S 4 ) 2 



V P=2 (T G ) = (1 - ^ 



+ 2S 2 (1 + S 4 )T G + 4S 4 T2 



(62) 



(63) 



((l-5 2 ) 2 + 45 2 T G ) 5 /2 

Our numerical results for a regular graph with two leads are plotted in figure 10 
together with the RMT results ( 6lj |62,|63). Notice that, the small conductances are 
emphasized, because of the presence of direct reflection and no direct transmission. At 
the same time, the most pronounced differences between the two symmetry classes are 
for Tq<1, where for (3 = 1 the conductance distribution diverges while it take a finite 
value for f3 = 2. 



6. Scattering from star-graphs. 

So far, we have studied scattering from well connected graphs, and we have shown 
that many statistical properties of the scattering matrix are described by RMT. We 
shall dedicate this section to demonstrate the statistical properties of the S— matrix for 
graphs which have non-uniform connectivity. 
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A representative example of this category are the "star" (or "Hydra") graphs 
14], [HJ. They consist of Vq bonds, all of which emanate from a single common vertex 
labeled with the index i = 0. The vertex at i = 0, will be referred to in the sequel as 
the head. The total number of vertices for such a graph is V = v + 1, and the vertices 
at the end of the bonds will be labeled by i — 1 • • • Vq. The star is a bipartite graph, a 
property which implies e.g., that there exists no periodic orbits of odd period [H|. To 
turn a star graph into a scattering system we add a lead to its head. 

It is a simple matter to derive the scattering matrix S = S^ M=1 ^ for a Neumann 
star. It reduces to a phase factor, 

The spectrum {k n } of the close system can be identified as the set of wave- vectors 
for which S(k) equals 1, which implies that no current flows in the lead. The resulting 
quantization condition is 

2£^ 1 tan(fc w Z, i 
££itaa(A;L<)+i 

which is satisfied once Yll=i tan(fc n Lj) = 0. This is identical with the condition derived 
in 



1 _ S (k) = <— > ^fZn r " = > ( 65 ) 



14 



The poles {n n } are the complex zeros of 

^^tan(K n Lj) + i = 0. (66) 
i=i 

To first order in — , we get 



r (1) = (67) 

£—/i=l cos(2fc n Li)+l 

which can be used as a starting point for the exact evaluation of the poles. For the 
latter one has to perform a self consistent search for the complex zeros of the secular 
equation (|66|). This is a time consuming process and the correct choice of the initial 
conditions is very important. 

In figure 11 we present our numerical results for the distribution of rescaled 
resonance widths V(j) for a star with = 20. The data are in excellent agreement 
with the RMT expectation given in Eq. (HI). We point that in this case the coupling to 
the continuum is weak since g ~ 10 ^> 1 and therefore the ^-distribution with M = 1 
is applicable. 

Using Eq. flBHj) we get the following relation for the scaled delay times 



A 2£ 



"0 



1 + (Z}i=i tan kLi) 

which can be used to generate V(T). The latter is reported in figure 12 together with 
the RMT prediction floTf). We notice that although the tail of the distribution agrees 
reasonably well with the RMT prediction, there are considerable deviations at the origin. 
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Figure 11. The rescaled resonance width distribution 7- > (7) for a star graph with 
v = 20. The solid line is the RMT prediction Eq. @). 



A peculiarity of the star graph is that the mean time delay is twice larger than 
the expected one from semiclassical considerations. To be more specific, for a generic 
chaotic system coupled to M channels, one has [EI] 



where (l p )k — T^ 1 is the average length of the classical paths inside the interaction area. 
Thus for the star with M — 1, we would expect based on Eq. (|69|) that (T)k — 2T c i. 
However, due to the fact that all the periodic orbits on the star graph are self tracing 
we get an additional factor of 2 and thus 

(T)* = l^d 1 - (70) 

The corresponding classical decay rate T c i can be found exactly by a direct 
evaluation of the eigenvalues of the classical evolution operator U. One can easily 
show that 




where flr(O) is the Vo x vq vertex scattering matrix at the star head as defined in Eq. (Eh, 
while 1 denotes the Vq x vq identity matrix. The square of the classical evolution operator 
U 2 has a block diagonal form 

U'-i^ 2 ° "| (72) 
and its spectrum consists of the values 1 — nq^p with multiplicity 2 and 1 — with 
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Figure 12. The distribution of the scaled partial delay times V(T) for a star graph 
with vq — 20. The dashed line is the RMT prediction Eq. (|58|). In the upper inset 
we show the same data in a double logarithmic plot. In the lower inset we plot the 
numerical results for < T >k (o). The solid line is the asymptotic value 1 expected 
from semiclassical considerations. 



multiplicity 2vq — 2. Therefore the spectrum of U is 



A » = ± i/ 1 - (TT%F (73) 



/ 4 

= ± Wl with multiplicity Vq — 1. 

\l l+v 

For short times where the classical evolution is applicable, the dominant eigenvalue is 



X u = \/l — leading to a classical decay rate 

2 



cl 



(74) 

l + ^o 

Substituting Eq. ([74]) we get eventually that (T% = v ^ob 1 which indicate that 
the mean time a particle spend inside the interaction regime is proportional to the 
Heisenberg time. Our numerical calculations reported in the lower inset of Fig. 12 
agrees nicely with the semiclassical prediction (|70|). 

Finally, we analyze the distribution of S when generated over different realizations 
of the lengths of the bond. For the one channel case this is equivalent with the 
distribution V{4>) of the phase of the 5-matrix. To derive the latter, it is convenient to 
rewrite the ^-matrix in the following form 

S = exp(i4>) = (-K + i)/(K + i), (75) 




Figure 13. The distribution of phases V(<f>) of the Si/-matrix for two stars with 
vq = 15 and vq = 30. The solid lines are the corresponding theoretical predictions 



with K = Y^l=i tan(fcXi). The probability distribution of K is 



V{K) = tan ( kL i) 



i=l 




AKx 



'dx 



AL 



^^g— ia; tan(fcL) 



Thus, with S = (S) 

V{4>) 



7TV 2 + K 2 

dK 



V(K) 



1-S 2 



(76) 



(77) 



Equation (|7 



2tt1 + 5 2 -25cos0 
reproduces the Poisson's kernel for a one-channel scattering matrix 



derived in the framework of RMT |07|. The conditions under which this result is derived 
in 



37 1 are fulfilled exactly in the present case. Our numerical results are reported in 
figure 13 and are in excellent agreement with the theoretical prediction Eq. (]77|) . 



7. Conclusions 

In this paper, we turned quantum graphs into scattering systems. We show that they 
combine the desirable features of both behaving "typically" and being mathematically 
simple. Thus, we propose them as a convenient tool to study generic behavior of chaotic 
scattering systems. 
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The classical dynamics on an open graph was defined, and the classical staying 
probability was shown to decay in an exponential way. The resulting classical escape 
rate was calculated and used to describe the properties of the corresponding quantum 
system. The scattering matrix was written in terms of classical orbits and an exact 
trace formula for the resonance density was found. A gap for the resonance widths 
has been obtained for "generic" graphs and its absence was explained for Neumann 
graphs. An analysis of the cross section autocorrelation function was performed and 
its non-universal characteristics were explained in terms of the short classical scattering 
trajectories. Finally, due to the relative ease by which a large number of numerical data 
can be computed for the graph models, we had performed a detail statistical analysis of 
delay times, resonance widths and distribution of the S'-matrix. Our results compares 
well with the predictions of RMT indicating that our model can be used in order to 
understand the origin of the connection between RMT and the underlying classical 
chaotic dynamics. 

The results reported here, complete our previous investigations on graphs. We 
conclude that quantum graphs may serve as a convenient paradigm in the area of 
quantum chaos, both for spectral and scattering studies. 
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